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Abstract. We present results of an investigation of clustering evolution of field galaxies be- 
tween a redshift of z ~ 1 and the present epoch. The current analysis relies on a sample of ~ 3600 
galaxies from the Calar Alto Deep Imaging Survey (CADIS). The multicolor classification and 
redshift determination is reliable up to / = 23 ma9 . The redshift distribution extends to z ~ 1.1, 
the resolution is Acz = 12000 km s . Thus the amplitude of the three-dimensional correlation 
function has to be estimated by means of the projected correlation function w(r p ). The validity 
of the deprojection was tested on the Las Campanas Redshift Survey (LCRS). The LCRS also 
serves as "local" measurement. We invented a new method to overcome the influence of redshift 
errors on w{rp). For evolution of the clustering strength the ansatz £(r com , z) oc (1 + z) q is used. 
For the galaxies as a whole the evolution parameter turns out to be q ~ —1.9, according to the 
prediction of linear theory. A formal dependency on the cosmology is presumably due to the 
small number of fields observed. However, the measured clustering growth clearly depends on 
Hubble type. At z ~ 1 early type galaxies are already much stronger clustered, an increase with 
q ~ — 1 is sufficient to explain the present day amplitude of the correlation function. 



1 The Calar Alto Deep Imaging Survey 

CADIS combines a very deep emission line survey carried out with an imaging Fabry-Perot inter- 
ferometer with a deep multicolour survey using three broad-band optical to NIR filters and up to 
thirteen medium-band filters when fully completed. The combination of different observing strategies 
facilitates not only the detection of emission line objects but also the derivation of photometric spectra 
of all objects in the fields without performing time consuming slit spectroscopy. 

For the multicolor sample, a classification scheme was developed, which is based on template spectral 
energy distributions (SEDs) ||. The classification algorithm basically compares the observed colors 
of each object with a color library of known objects. This color library is assembled from observed 
spectra by synthetic photometry performed on an accurate representation of the instrumental charac- 
teristics used by CADIS. Using the minimum variance estimator, each object is assigned a type (star 
- QSO - galaxy), a redshift (if it is not classified as star), and an SED (in terms of a number, = 
E0, 100 = starburst). The formal errors in this process depend on magnitude and type of the object 
and are of the order of a z = 0.017, and cjsed = 2, respectively. 

The seven CADIS fields measure ~ 1/30 D° each and are located at high Galactic latitude to avoid 
dust absorption and reddening. Four of these have been fully analysed so far. We identified 3626 
galaxies with I < 23 in the redshift range 0.2 < z < 1.07. 



2 The projected correlation function 

The projected correlation function w(r p ) was first introduced by to overcome the influence of non- 
negligible peculiar velocities on the three-dimensional correlation function. Redshift inaccuracies also 
increase the noise and supress the correlation signal. Therefore we used the projected correlation 
function to investigate the evolution of galaxy clustering. It is defined by 
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r p is the projected distance between pairs of galaxies (the distance perpendicular to the line of sight), ir 
is the distance between the two galaxies parallel to the line of sight, and r 2 = r 2 -|-7r 2 . If £(r) = (?"/Vo)~ 7 , 
then equation @) yields 

w(r p ) = C r y p ^ (2) 

with 
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Thus computing w(r p ) provides a measurement of the parameters of the three-dimensional correlation 
function, namely ro and 7. 

Following [|J one can calculate the projected correlation function from 
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Since the three-dimensional two-point correlation function has the form of a power law, it converges 
rapidly to zero with increasing pair separation. Therefore the integration limits do not have to be 
±00, they only have to be large enough to include all correlated pairs. Since the observable in the 
first place is the redshift z instead of the physical separation, we make a coordinate transform [B]: 
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for ft A = 0. 

The way to estimate w(r p ) in practice is to count the projected distances between pairs of galaxies 
that are separated in redshift space by not more than 5z, in appropriate projected distance bins. We 
use the estimator by |§, which we will call Cesti(f p ) in the following. To derive w(r p ), Cesti(?>) has to 
be multiplied with the "effective depth" Ary in which galaxies are taken into account. 
In reality, one has to cope with a selection function of some kind or another, and not with a top-hat 
function (of probability unity to find a galaxy within the borders of the survey and zero otherwise). 
The varying probability to find pairs of galaxies separated by a redshift 5z within the survey has to 
be taken into account. It can be included in the calculation by multiplying the integrand in equation 
(||) with the squared redshift distribution, normalised to unity at it's maximum. With this correction 
for the selection function, the "depth" converges to a fixed value for 5cz —* 00 and does not grow 
anymore, even if the integration limits cover more than the total depth of the survey, thus 
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for n A = 0. 

The choice of the integration limit 5 Z is somewhat arbitrary. To find the appropriate integration lim- 
its, we simulated the influence of redshift errors on the projected correlation function. We added 
artificial redshift errors to redshifts of LCRS galaxies. The errors (one time for a resolution of 
Acz = 5000 km s" 1 and the other time with Acz = 12000 km s _1 , the size of the CADIS redshift 
resolution) were randomly drawn out of a Gaussian error distribution. We calculated the projected 
correlation function for the modified samples for increasing integration limits 8cz. The amplitude was 
fitted at r p = 500/i _1 kpc, in the range —1.15 < logr p < —0.3, to make sure we fit it where the signal 
to noise is high. Also the choice of r p = 500/i _1 kpc facilitates the comparison with the CADIS data, 
where we fit the amplitude at r p ~ 316/i _1 kpc (at the mean redshift of the survey (z) this corresponds 
to a comoving separation of ~ 505/i -1 kpc). 

Figure [l] shows the amplitude of w(r p = 500/i~ 1 kpc), for increasing integration limits. The measure- 
ment is compared to the modified data. 
One starts to sample the correlation signal when the integration limits are larger than the full width 
at half maximum of the redshift error distribution. Also the amplitude is diminished. For errors of 
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Figure 1: The influence of redshift measurement errors on the projected two-point correlation function, 
shown is the weighted mean of four sectors. w{r p ) for increasing integration limits is calculated for 
an Acz = 5000 km s _1 (dotted line), and for Acz = 12000 km s _1 (dashed line). The errors of the fits 
are not plotted here to avoid confusion. The marks indicate the integration limits which have to be 
chosen for the calculation of w{r p ). 

Acz = 12000 km s -1 3 the maximum amplitude is a factor 1.4 lower than in the case of the unchanged 
data. 

We calculated w(r p ) for the CADIS data with ±5cz = 15000 km s . To facilitate the direct comparison 
with the LCRS data, we used the modified sample (with artificial errors of size Acz = 12000 km s _1 ) 
and the same integration limits. 

Bright stars in our field were masked out, and a random catalogue consisting of 10000 "galaxies" was 
generated with the same properties as the real data. The calculation was carried out for a closed 
high-density model (Qq = 1, ft\ = 0), a hyperbolic low-density model (f2o = 0.2, S7a = 0), and a 
flat low-density model with non-zero comological constant (Oq = 0.3, ft\ = 0.7). Figure [2] shows the 
results for the latter two ones. 
From the amplitude of w(r p ) we can deduce the amplitude of the three-dimensional correlation func- 
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Figure 2: Projected correlation function in three redshift bins, for an open model (left panel), and for 
a flat universe with a non-zero cosmological constant (right panel). 
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with the numerical factor C from equation (||). 

We can now parametrise the evolution with a parameter q which gives directly the deviation from the 
global Hubble flow: 

£(r com = IMpc) = &(!+*)* . (8) 

q can be deduced straightforwardly: if one fits log£(r com = IMpc) versus log(l + z), it is just the slope 
of the straight line. 

We find q = -2.68 ± 0.16 for £l = 1, n A = 0, q = -1.92 ± 0.17 for tt = 0.2, Q A = 0, and 
q = —1.23 ± 0.20 for 0,q = 0.3, £l\ = 0.7. The formal dependence on the cosmological model adopted 
for the calculation is larger than expected from a simple examination of the distances corresponding 
to a certain angle 0. Therefore we regard this as an indication that the differences are presumably 
due to the small number of fields observed. 

We can compare our results for the Qq = 1, S7a = case with the results of ||. For the direct 
comparison we have to multiply our measured amplitudes of the projected correlation function by 1.4 
to correct for the influence of large redshift errors. 

With this correction, the CFRS data points are consistent with our own measurement, although with 
large errors, see Figure ||. Q claim that if ro(z = 0) = 5/i _1 Mpc, < e 2. The fit of their data 
points, including the connection to z = 0, yields q = —3.043 ± 0.213 (= e = 1.8). If the connection 
to z = is disregarded, we find q = —1.184 ± 0.634. The fit including our redshift error corrected 
LCRS point instead of the 0] point yields q = —2.298 ± 0.238. This is even a bit less than our own 
measurement (q = —2.68 ± 0.16), but nevertheless equal within the errors. 

2.1 The evolution of clustering for different Hubble types 

Our multicolor classification scheme gives the SED of each galaxy, and thus we are able to calculate 
the projected correlation function for different Hubble types. We divided our sample at SED = 60, 
the sample with < SED < 60 is called "early type", and the sample with 60 < SED < 100 is called 
"late type" in the following. The amplitudes of the correlation function of the early type galaxies 
are at all redshifts significantly larger than for the late type ones, and also clearly larger than for the 
complete CADIS sample. Also the evolution of the clustering strength is different. The amplitude 
of the correlation function of the early type galaxies evolves much slower, we find q = —1.71 ± 0.39 
for Q = 1, Q A = 0, q = -1.37 ± 0.36 for VL Q = 0.2, ft A = 0, and q = -0.67 ± 0.33 for Q = 0.3, 
f^A = 0.7. For the late type sample the errors are too large to facilitate a quantitative statement about 
the custering evolution. 

Combining SED and redshift information, we are able to calculate restframe B band luminosities 
for all the galaxies. Therefore we divided our sample at Mb = —I8 ma9 , and calculated the projected 
correlation function for bright and faint galaxies in the redshift range 0.3 < z < 0.6, where the number 
of faint galaxies is sufficient to enable the determination of the amplitude of the projected correlation 
function. The bright galaxies are clearly higher clustered than the faint ones. Figure ^ shows the 
amplitudes of the three-dimensional correlation function and the fit for the evolution paramter q for 
the complete sample, for early and late type galaxies and also for the bright and the faint sample. 
The left panel shows the results for f?o = 0.3, S7a = 0.7, the right panel shows the comparison of our 
data with the CFRS data. 



3 Conclusion 

We find that the amplitude of the correlation function of all the galaxies in our sample grows with 
q ~ —1.9 between a redshift of z ~ 1.1 and the present epoch. The rate of clustering growth is 
consistent with the results of linear perturbation theory {q = —2). The amplitude of the early type 
galaxies grows much slower with redshift (q ~ —1)- This result can be explained in the context of 
biased galaxy formation, if we assume a substantial evolution of the galaxies between a redshift of 
z ~ 1.1 and today. The first generation of galaxies forms in a highly clustered state, in the bumps and 
wiggles superimposed on the very large scale perturbations of the dark matter density field, whereas 
the next generations of galaxies form later in the lower clustered environment in the wings of the super 
large scale overdensities. Therefore early type galaxies are much stronger clustered than the young, 
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Figure 3: The evolution of the clustering strength (at lh~ l Mpc) with redshift. Left panel: The line 
is the fit to the data points, the first data point is the weighted mean of four LCRS sectors, the three 
other ones are CADIS data, for a flat universe with a non-zero cosmological constant. The LCRS 
point is not included in the fit for the early type galaxies. Right panel: The amplitudes of the three 
dimensional correlation function at r com = l/i -1 Mpc, deduced from the projected correlation function 
of the CFRS data (crosses), in comparison with our own data (filled symbols). The open symbols are 
our data points corrected for the influence of the redshift errors on the amplitude of the projected 
correlation function. The dotted line is the fit of the CFRS data points including the value of tq(z = 0) 
from 0, the dashed line is the fit using the corrected LCRS point instead. 



late type ones even at higher redshifts. While the universe evolves and the clustering of the underlying 
dark matter density field grows, galaxies age, and eventually merge to form massive ellipticals. They 
add to the early type population at different times, thus the early type sample consists of galaxies 
of different ages, which have formed in increasingly less clustered states. This means that biasing 
decreases with increasing redshift. The net effect is a very slow growth of clustering of early type 
galaxies. 

Merging also plays an important role, because of two effects on the correlation function: first of all, 
galaxies change their SEDs on relatively small timescales. They "suddenly" disappear from the lower 
clustered late type sample and reappear in the early type sample. The second effect is that while two 
galaxies merger to form one, the small distances disappears, the probability of finding a galaxy near 
another galaxies decreases, and therefore the amplitude of the correlation function decreases. 
If bright galaxies form in the highest peaks of the dark matter field, they are expected to be higher 
clustered than the faint ones. The difference in clustering strength is larger for the early type/late type 
sample than for bright and faint galaxies, which corroborates the hypothesis, that galaxy evolution 
and merging plays an important role in the evolution of clustering. 
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